Observation of B — > ppK 

K. Abe 9 , K. Abe 40 , R. Abe 30 , T. Abe 41 , I. Adachi 9 , Byoung Sup Ahn 17 , H. Aihara 42 , M. Akatsu 23 , Y. Asano 47 , 
T. Aso 46 , V. Aulchenko 2 , T. Aushev 14 , A. M. Bakich 37 , Y. Ban 34 , E. Banas 28 , S. Behari 9 , P. K. Behera 48 , 
A. Bondar 2 , A. Bozek 28 , M. Bracko 21 ' 15 , J. Brodzicka 28 , T. E. Browder 8 , B. C. K. Casey 8 , P. Chang 27 , Y. Chao 27 , 
K.-F. Chen 27 , B. G. Cheon 36 , R. Chistov 14 , S.-K. Choi 7 , Y. Choi 36 , M. Danilov 14 , L. Y. Dong 12 , J. Dragic 22 , 
A. Drutskoy 14 , S. Eidelman 2 , V. Eiges 14 , Y. Enari 23 , C. W. Everton 22 , F. Fang 8 , H. Fujii 9 , C. Fukunaga 44 , 
M. Fukushima 11 , N. Gabyshev 9 , A. Garmash 2 ' 9 , T. Gershon 9 , K. Gotow 49 , R. Guo 25 , J. Haba 9 , H. Hamasaki 9 , 
F. Handa 41 , K. Hara 32 , T. Hara 32 , H. Hayashii 24 , M. Hazumi 9 , E. M. Heenan 22 , I. Higuchi 41 , T. Higuchi 42 , 
T. Hojo 32 , T. Hokuue 23 , Y. Hoshi 40 , S. R. Hou 27 , W.-S. Hou 27 , H.-C. Huang 27 , T. Igaki 23 , Y. Igarashi 9 , H. Ikeda 9 , 
K. Inami 23 , A. Ishikawa 23 , H. Ishino 43 , R. Itoh 9 , H. Iwasaki 9 , Y. Iwasaki 9 , P. Jalocha 28 , H. K. Jang 35 , J. H. Kang 50 , 
J. S. Kang 17 , P. Kapusta 28 , N. Katayama 9 , H. Kawai 3 , H. Kawai 42 , Y. Kawakami 23 , N. Kawamura 1 , T. Kawasaki 30 , 
H. Kichimi 9 , D. W. Kim 36 , H. J. Kim 50 , H. 0. Kim 36 , Hyunwoo Kim 17 , S. K. Kim 35 , K. Kinoshita 5 , H. Konishi 45 , 
S. Korpar 21 * 15 , P. Krizan 20 ' 15 , P. Krokovny 2 , S. Kumar 33 , A. Kuzmin 2 , Y.-J. Kwon 50 , J. S. Lange 6 , G. Leder 13 , 

S. H. Lee 35 , A. Limosani 22 , D. Liventsev 14 , R.-S. Lu 27 , J. MacNaughton 13 , G. Majumder 38 , F. Mandl 13 , 
T. Matsuishi 23 , S. Matsumoto 4 , Y. Mikami 41 , W. Mitaroff 13 , K. Miyabayashi 24 , Y. Miyabayashi 23 , H. Miyake 32 , 
H. Miyata 30 , G. R. Moloney 22 , T. Nagamine 41 , Y. Nagasaka 10 , Y. Nagashima 32 , M. Nakao 9 , J. W. Nam 36 , 
Z. Natkaniec 28 , S. Nishida 18 , O. Nitoh 45 , S. Noguchi 24 , T. Nozaki 9 , S. Ogawa 39 , T. Ohshima 23 , T. Okabe 23 , 
S. Okuno 16 , S. L. Olsen 8 , W. Ostrowicz 28 , H. Ozaki 9 , P. Pakhlov 14 , H. Palka 28 , C. S. Park 35 , C. W. Park 17 , 
K. S. Park 36 , L. S. Peak 37 , J.-P. Perroud 19 , M. Peters 8 , L. E. Piilonen 49 , F. Ronga 19 , N. Root 2 , K. Rybicki 28 , 

H. Sagawa 9 , Y. Sakai 9 , H. Sakamoto 18 , M. Satapathy 48 , A. Satpathy 9 ' 5 , O. Schneider 19 , S. Schrenk 5 , 
C. Schwanda 9 ' 13 , S. Semenov 14 , K. Senyo 23 , H. Shibuya 39 , B. Shwartz 2 , V. Sidorov 2 , J. B. Singh 33 , S. Stanic 47 , 
A. Sugi 23 , A. Sugiyama 23 , K. Sumisawa 9 , K. Suzuki 9 , S. Y. Suzuki 9 , T. Takahashi 31 , F. Takasaki 9 , M. Takita 32 , 
K. Tamai 9 , N. Tamura 30 , J. Tanaka 42 , G. N. Taylor 22 , Y. Teramoto 31 , S. Tokuda 23 , T. Tomura 42 , S. N. Tovey 22 , 
K. Trabelsi 8 , T. Tsuboyama 9 , T. Tsukamoto 9 , S. Uehara 9 , K. Ueno 27 , Y. Unno 3 , S. Uno 9 , K. E. Varvell 37 , 
C. C. Wang 27 , C. H. Wang 26 , J. G. Wang 49 , M.-Z. Wang 27 , Y. Watanabe 43 , E. Won 35 , B. D. Yabsley 9 , Y. Yamada 9 , 
M. Yamaga 41 , A. Yamaguchi 41 , H. Yamamoto 41 , Y. Yamashita 29 , M. Yamauchi 9 , S. Yanaka 43 , P. Ych 27 , 
M. Yokoyama 42 , Y. Yuan 12 , J. Zhang 47 , Y. Zheng 8 , and D. Zontar 47 



(Belle Collaboration) 

1 Aomori University, Aomori 
2 Budker Institute of Nuclear Physics, Novosibirsk 
3 Chiba University, Chiba 
4 Chuo University, Tokyo 
5 University of Cincinnati, Cincinnati OH 

6 University of Frankfurt, Frankfurt 
7 Gyeongsang National University, Chinju 
8 University of Hawaii, Honolulu HI 
9 High Energy Accelerator Research Organization (KEK), Tsukuba 
10 Hiroshima Institute of Technology, Hiroshima 
11 Institute for Cosmic Ray Research, University of Tokyo, Tokyo 
12 Institute of High Energy Physics, Chinese Academy of Sciences, Beijing 
13 Institute of High Energy Physics, Vienna 
14 Institute for Theoretical and Experimental Physics, Moscow 
15 J. Stefan Institute, Ljubljana 
16 Kanagawa University, Yokohama 
17 Korea University, Seoul 
18 Kyoto University, Kyoto 
19 IPHE, University of Lausanne, Lausanne 
20 University of Ljubljana, Ljubljana 

21 University of Maribor, Maribor 
22 University of Melbourne, Victoria 
23 Nagoya University, Nagoya 
24 Nara Women's University, Nara 

25 National Kaohsiung Normal University, Kaohsiung 

26 National Lien-Ho Institute of Technology, Miao Li 

27 National Taiwan University, Taipei 
28 H. Niewodniczanski Institute of Nuclear Physics, Krakow 



1 



29 Nihon Dental College, Niigata 

30 Niigata University, Niigata 
31 Osaka City University, Osaka 

32 Osaka University, Osaka 
33 Panjab University, Chandigarh 
34 Peking University, Beijing 
Seoul National University, Seoul 
36 Sungkyunkwan University, Suwon 
37 University of Sydney, Sydney NSW 
38 Tata Institute of Fundamental Research, Bombay 
39 Toho University, Funabashi 
40 Tohoku Gakuin University, Tagajo 
41 Tohoku University, Sendai 
42 University of Tokyo, Tokyo 
43 Tokyo Institute of Technology, Tokyo 
44 Tokyo Metropolitan University, Tokyo 
45 Tokyo University of Agriculture and Technology, Tokyo 
46 Toyama National College of Maritime Technology, Toyama 
47 University of Tsukuba, Tsukuba 
48 Utkal University, Bhubaneswer 
49 Virginia Polytechnic Institute and State University, Blacksburg VA 
50 Yonsei University, Seoul 



We report the observation of the decay mode B — ► ppK based on an analysis of 29.4 fb~ 
of data collected by the Belle detector at KEKB. This is the first example of a b — > a transition 
with baryons in the final state. The pp mass spectrum in this decay is inconsistent with phase 
space and is peaked at low mass. The branching fraction for this decay is measured to be B(B ± — + 
ppK^ = (4.31^9 (stat) ^ 0.5(syst)) x 1CP 6 . We also report upper limits for the decays B° — > ppK$ 
and B ± — > ppn ± . 



PACS numbers: 13.20.He, 13.25.Hw, 13.60.Rj 
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We report the results of searches for the decay modes 
B + — ► ppK + H and B° — > ppKg. These modes are 
expected to proceed mainly via b s penguin diagrams. 
We also search for B + — ► ppn + which is expected to 
occur primarily via a b — > u tree process. Once they are 
established, these baryonic modes may be used to either 
constrain or observe direct CP violation in B decay 

In contrast to charm meson decay, final states with 
baryons are allowed in B meson decay. To date, a few 
low multiplicity B decay modes with baryons in the final 
state from b — ► c transitions have been observed || . Rare 
B decays due to charmless b —> s and b — > u transitions 
should also lead to final states with baryons. A num- 
ber of searches for such modes have been carried out by 
CLEO j§, ARGUS §, and LEP § but only upper lim- 
its were obtained. Stringent upper limits for two-body 
modes such as B° — + pp, B + — ► Ap and B° — ► AA have 
recently been reported by Belle . 

We use a 29.4 fb _1 data sample, which contains 31.9 
million produced BB pairs, collected with the Belle de- 
tector at the KEKB asymmetric-energy e + e~ (3.5 on 
8 GeV) collider ||. KEKB operates at the T(45) res- 
onance (y/s — 10.58 GeV) with a peak luminosity that 
exceeds 5 x 10 33 cm~ 2 s . The Belle detector is a 
large-solid-angle magnetic spectrometer that consists of 
a three-layer silicon vertex detector (SVD), a 50-layer 
central drift chamber (CDC), a mosaic of aerogel thresh- 
old Cerenkov counters (ACC), time-of- flight scintillation 
counters (TOF), and an array of CsI(Tl) crystals (ECL) 
located inside a superconducting solenoid coil that pro- 
vides a 1.5 T magnetic field. An iron flux-return located 
outside of the coil is instrumented to identify Kl and 
muons (KLM). The detector is described in detail else- 
where ||. 

We select well measured charged tracks with impact 
parameters with respect to the interaction point of less 
than 0.3 cm in the radial direction and less than 3 cm 
in the beam direction (z). These tracks are required to 
have p T > 50 MeV/c. 

Particle identification likelihoods for each particle hy- 
pothesis are calculated by combining information from 
the TOF, ACC system with dE/dx measurements in the 
CDC. Protons and anti-protons are identified using all 
particle ID systems and are required to have proton like- 
lihood ratios (L P /(L P + Lk) and L P /(L P + L^)) greater 
than 0.6. Proton candidates that are electron-like accord- 
ing to the information recorded by the CsI(Tl) calorime- 
ter are vetoed. This selection is 89% efficient for protons 
with a 7% kaon misidentification rate. To identify kaons 
(pions) , we require the kaon (pion) likelihood ratio to be 
greater than 0.6. This requirement is 88% efficient for 
kaons with a 8.5% misidentification rate for pions. In 
addition, we remove kaon candidates that are consistent 
with being protons. 



For the B — > ppK$ mode, we select Kg candidates 
from tt + tt~ candidates that lie within the mass window 
0.482 GeV/c 2 < M(7T+7r-) < 0.514 GeV/c 2 (±4cr). The 
distance of closest approach between the two daughter 
tracks is required to be less than 2.4 cm. The impact pa- 
rameter of each track in the radial direction should have 
magnitude greater than 0.02 cm, and the flight length 
should be greater than 0.22 cm. The difference in the 
angle between the vertex direction and the K$ flight di- 
rection in the x — y plane is required to satisfy Atp < 0.03 
rad. 

To reconstruct signal candidates in the B + — > ppK + 
mode, we form combinations of a kaon, proton and anti- 
proton that are inconsistent with the following b — > ccs 
transitions: B + — * J/i/jK + , J/tp — > pp; B + — > i] c K + , 
rj c -> pp: B+ -> tp'K + , ip' -> pp and B+ -> Xc[o,i] K+ , 
Xc[o,i] — > PP- This set of requirements is referred to as 
the charm veto ||To|| . Similar charm vetoes are applied 
in the analysis of the other decay modes. In the case of 
B° — > ppKs, events with pK$ or pK$ masses consistent 
with the A c are rejected. 

To isolate the signal, we form the beam-constrained 

mass, Mbc = \Awm - Brecon! and energy difference 
AE = -Brecon — -E-beam hi the T(45') center of mass frame. 
Here E hclim , E rccon and P locon are the beam energy, the 
reconstructed energy and the reconstructed momentum 
of the signal candidate, respectively. The signal region 
for AE is ±50 MeV which corresponds to ±5cr where a 
is the resolution determined from a Gaussian fit to the 
Monte Carlo (MC) simulation. The signal region for M DC 
is 5.270 GeV/c 2 < M hc < 5.290 GeV/c 2 . The resolution 
in beam-constrained mass is 2.8 MeV/c 2 and is domi- 
nated by the beam energy spread of KEKB. 

Several event topology variables provide discrimina- 
tion between the large continuum (e + e~ — > qq, where 
q = u,d,s,c) background, which tends to be collimated 
along the original quark direction, and more spherical 
BB events. We form a likelihood ratio using two vari- 
ables. Six modified Fox- Wolfram moments and the cosine 
of the thrust angle are combined into a Fisher discrimi- 
nant pi]] . For signal MC and continuum data, we then 
form probability density functions for this Fisher discrim- 
inant, and the cosine of the B decay angle with respect to 
the z axis (cos 0b)- The signal (background) probability 
density functions are multiplied together to form a signal 
(background) likelihood Cs (Cbg)- The likelihood ratio 
£s/(£s + £bg) is then required to be greater than 0.6. 
The event topology requirements retain 78% of the signal 
while removing 87% of the continuum background. 

In Fig. 0, we show the AE and beam-constrained mass 
distributions for the signal candidates. We fit the AE dis- 
tribution with a double Gaussian for signal and a linear 
background function with slope determined from the My, c 
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sideband. The mean of the Gaussian is determined from 
B" — > A c p7r+7r~, A c — > pK^n + decays. The fit to the 
AE distribution gives a yield of 42.8±g° 6 8 with a signifi- 
cance of 5.6 standard deviations. In the fit to the AE dis- 
tribution, the region with AE < — 120 MeV is excluded 
to avoid feed-downs from modes such as B — ► ppK* . As 
a consistency check, we fit the M\, c distribution to the 
sum of a signal Gaussian and a background function with 
kinematic threshold. The width of the Gaussian is fixed 
from MC simulation while the mean is determined from 
B + — ► D°tt + data. The shape parameter of the back- 
ground function is determined from AE sideband data. 
In the Mbc distribution, we observe a signal of 42.9^9 "f 
events. The signal yields and the branching fractions are 
determined from fits to the AE distribution rather than 
Afbc to minimize possible biases from BB background 
which tends to peak in Mx> c but not in AE. 




FIG. 1. (a) AE and (b) M bc distributions for B + -> ppK 



candidates. 

The background in these modes is predominantly due 
to continuum events. To check for BB backgrounds 
that might peak in the signal region, we used two large 
BB MC samples that correspond to an integrated lu- 
minosity that is about twice the size of the data sam- 
ple. The estimated background is of order one event 
and no backgrounds that peak in the AE signal region 
were found. We also examined dedicated MC samples of 
b — ► c decay modes with baryons in the final state. We 
restricted our attention to low multiplicity decay modes. 
We generated samples of B° — > A+p, B~ — ► A+p7r~ and 
B~ — ► k+pe~v e that correspond to an integrated lumi- 
nosity about a factor of ten larger than the data sample 
used here. The A c charmed baryon was allowed to de- 
cay into all measured decay modes that contain a proton. 
Again no peaking backgrounds were observed. 

We also examine the M{pp) mass distributions for 
events in the AE, Mb c signal region. The signal yield 
as a function of pp mass is shown in Fig. |[ These yields 
were determined by fits to the AE distribution in bins of 
pp invariant mass. The distribution from a three-body 
phase space MC normalized to the area of the signal is 
superimposed. It is clear that the observed mass distri- 



bution is not consistent with three-body phase space but 
instead is peaked at low pp mass. We also examine the 
pK~ mass distribution but do not observe any obvious 
narrow structures such as the A(1520). 

To avoid model dependence in the determination of the 
branching fraction for ppK + , we fit the AE signal yield 
in bins of M(pp) and correct for the detection efficiency 
in each bin using a three-body phase space B + — ► ppK + 
MC model. The results of the fits are given in Table I. 
We then sum the partial branching fractions in each bin 
to obtain 

B{B+ -► ppK+) = (4.3tg g(stat) ± 0.5(syst)) x 1CT 6 . 

For M(pp) < 3.4 GeV/c 2 , the mass region below the % c 
and ip' resonances, B(B+ -> ppK+) = (4.4tj^(stat) ± 
0.5(syst)) x 10~ 6 with the charm veto applied. For 
M(pp) < 2.8 GeV/c 2 , the region below charm threshold, 
we obtain B{B+ -> ppK+) = (3.9±^(stat)±0.4(syst)) x 
1(T 6 . 



TABLE I. Fit results in bins of M(pp). The detection effi- 
ciency (tdetect) and the partial branching fraction (£>) for each 
bin are also listed. 



M(pp)(GcV/c 2 ) 


AE yield 


^detect 


£>(x 10 -6 ) 


< 2.0 
2.0-2.2 
2.2-2.4 
2.4-2.6 
2.6-2.8 
2.8-3.4 
3.4-4.0 
4.0-4.8 


10.2+1* 

11 q+4.6 
rx.»_ 3 9 

c c+3.7 
d - J -3.0 

3.3l2;J 

1 n+2.5 
1 - z '-2.2 

^+ 3 - 5 

U.O_ 2 g 


0.33 
0.34 
0.30 
0.29 
0.30 
0.29 
0.27 
0.25 


Q7+ tJ - 4y 

u - y '-0.35 
n 7 n+0.39 
u - ' J -0.32 

1 24+ ' 48 
61 + ' 41 

u - Di -0.33 

34+0-32 

o ™+ 38 

U.OU_ 29 

04+ ' 45 



The contributions to the systematic error for the B + — > 
ppK + mode are the uncertainties due to the tracking ef- 
ficiency (6%), particle identification efficiency (8%) and 
the modeling of the likelihood ratio cut (2.6%). The 
particle identification systematic includes contributions 
of 3% for the proton and anti-proton and 2% for the 
charged kaon. The error in proton/anti-proton identifi- 
cation is determined using A/A samples, while the error 
in kaon identification efficiency is obtained from kinemat- 
ically selected D* + -> D°ir + ,D° -> R-tt+ in the data. 
The systematic error due to the modeling of the likeli- 
hood ratio cut is determined using B+ -> D°tt+ events 
reconstructed in data. The systematic error in the yield 
of the AE fit (3.8%) was determined by varying the mean 
and a of the signal and the shape parameter of the back- 
ground. The sources of systematic error are combined in 
quadrature to obtain the final systematic error of 11.0%. 

For events in the AE, M^ c signal region we examine 
the proton, anti-proton and kaon particle identification 
likelihood distributions and compare to signal MC simu- 
lation. No discrepancy is observed. We also verify that 
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FIG. 2. The fitted yield divided by the bin size for 
B + — > ppK + as a function of pp mass. The charm veto is 
applied. The distribution from non-resonant B + — > ppK + 
MC simulation is superimposed. The inset shows the pp mass 
distribution for the J/ipK + signal region. 



the ECL shower width distribution is consistent with MC 
expectations for the proton and anti-proton candidates. 
In addition, we check the branching fraction as the cuts 
on the proton and anti-proton probabilities and likeli- 
hood ratio are varied. We do not observe any systematic 
trends beyond statistics. 

To verify the analysis procedure and branching frac- 
tion determination, we remove the J/ip veto and ex- 
amine the decay chain B + — > J/i/jK + , J/ip — > pp. A 
clear signal of 26.4 ± 5.2 events is then observed in the 
AE spectrum. We also observe 25.9 ± 5.1 events in 
the Mbc distribution. The pp invariant mass spectrum 
for JJjl)K + signal candidates is shown as an inset in 
Fig. g. We use the AE yield and the MC detection 
efficiency of 0.30 to determine the branching fraction 
B(B+ -► J/iPK+) = (13.1 ± 2.6) x 10~ 4 . This is in 
good agreement with the PDG world average, B(B + — > 
J/ipK+) = (10.0±1.0)xl0- 4 @, which was obtained by 
experiments that reconstruct the J/tp in dilepton modes. 

We also examined two related decay modes B° — > 
ppKs and B + — ► ppir^ that may help clarify the inter- 
pretation of the signal. Measurement of B° — > ppKs will 
help to determine the role of the spectator quark in b — > s 
decays with baryons, while observation of B + — > ppn + 
will constrain the ratio of the b — > u tree and b — > s 
penguin diagrams in decays with baryons. 

For B° — > ppKs, after the application of the charm 



events. Applying the Feldman- 
we obtain an upper limit of 



and A c vetoes, no significant signal is observed in ci- 
ther the AE or Mb c distribution. A fit to the AE dis- 
tribution gives 6.4+ 4 ; 4 - 
Cousins procedure |13fl 
less than 16 events at the 90% confidence level (C.L.). 
After reducing the detection efficiency by the system- 
atic error, we obtain an upper limit at 90% C.L. of 
B(B° ^ppK°) < 7.2 x 10- 6 . 

In the B + — > ppit + mode, after the application of 
the charm veto we perform a fit to the AE distribution 
that allows for B + — > ppir + signal and a reflection from 
misidentified B + — > ppK + decays. This fit gives a signal 
yield of 16.2 jl|'Q events and a significance of 2.1c. The 
excess in the AE fit corresponds to a branching fraction 
B(B+ -> ppir+) = (l-9±o - 9±0.3) x 10~ 6 or an upper limit 
of B{B+ -> ppn+) < 3.7 x 10~ 6 at 90% C.L. after taking 
into account the systematic error. 

We have observed a significant signal (5.6<j) for the 
decay B + — * ppK + . This is the first b — * s decay mode 
with baryons in the final state. In the future, this mode 
can be used to search for direct CP violation j2|. Wc 
find that its pp mass spectrum is inconsistent with phase 
space and is peaked toward low mass. This feature is sug- 
gestive of quasi two-body decay. It is also possible that 
the decay is a genuine three-body process and that this 
feature of the M{pp) spectrum is a baryon form factor 
effect @|n§. 
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